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Abstract. We present models of the spectral reflectances of groups of outer Solar 
System objects defined primarily by their colors in the spectral region 0.4 -1.2 (im , 
and which have geometric albedo 0.04 at wavelength 0.55 fim . Our models of the 
groups with the strongest reflectance gradients (reddest colors) use combinations of 
organic tholins. We test the hypothesis that metal-reddened igneous rock-forming 
minerals contribute to the red colors of Centaurs and KBOs by using the space- 
weathered lunar soil as one of the components of our models. We find that our 
models can admit the presence of moderate amounts of space-weathered (metal- 
reddened) minerals, but that they do not require this material to achieve the red 
colors of the reddest outer Solar System bodies. Our models with organic tholins 
axe consistent with the results of other investigators. 


1. INTRODUCTION 


The compositions of the small bodies in the outer Solar System (OSS), 
particularly those with geometric albedo = 0.1 at A = 0.5 /xm , are un- 
derstood only in broad outline. While many outer belt asteroids, small 
planetary satellites, Centaur objects, Kuiper Belt objects (KBOs), and 
comet nuclei have been observed photometrically and spectroscopically 
(primarily at wavelengths 0.3 < A < 2.5 /jm ),. relatively little direct 
information on their compositions has resulted. The natural materials 
expected to compose these bodies are ices of H 2 O and other volatile 
molecules, minerals, and solid macromolecular carbonaceous material. 
The term ’’macromolecular carbonaceous” includes, but is not limited 
to, the complex organic materials found in the carbonaceous mete- 
orites (Cronin etal. 1988, Cronin and Chang 1993), interplanetary 
dust particles (Bradley etal. 1996), and micrometeorites (Maurette 
etal. 1996). It also includes complex organic materials presumed to 
have been synthesized by the deposition of energy in CHON-rich matter 
directly on the surfaces of KBOs and other solid bodies in the outer 
Solar System. 

Absorption bands of H 2 O ice have been detected in a few Centaurs 
(e.g., Brown et al. 1998), Kuiper Belt objects (e.g., Brown etal. 1999), 
irregular planetary satellites (Owen et al. 1999), and at least one comet 
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Wavelength (pm ) 

Figure 1. Four groups defined by photometry of Centaurs and Kuiper Belt objects 
by Barucci etal. (2001). Representative members of the groups axe: Group I (1996 
T066, 2060 Chiron), Group II (1997 CU26, 1998 SG35), Group III (1993 SC, 1997 
CQ29), Group IV (5145 Pholus, 1994 TB). We have normalized these curves to 
geometric albedo 0.04 at wavelength 0.56 /rm . 

3. MODEL COMPONENTS 
3.1. Organic Solids 

We are especially concerned in this paper with Solar System bodies 
having large reflectivity gradients, or red colors. When the very red 
Centaur object 5145 Pholus was first observed in the 0.3-1.0^m region, 
three different investigator teams (Fink et al. 1992; Binzel 1992; Davies 
et al. 1993) all concluded that only solid organic materials, or materials 
similar to a solid polymer of hydrogen cyanide, could account for the red 
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color. Simple spectrum matching with the reflectance of the synthetic 
organic material called tholin provided a reasonably good fit to the color 
of Pholus in that limited wavelength region, while plausible minerals 
and ices did not. 

The Trojan asteroid 624 Hektor (classified as D) is among the red- 
dest of the traditional asteroids, although it is much less red than the 
reddest of the KBOs and Centaurs. Cruikshank etal. (2001) showed 
that the spectral reflectance of Hektor can be modeled with Mg-rich 
pyroxene and serpentine, over the wavelength range 0.3 - 4 jim . Their 
result demonstrated that organic solids are not essential for modeling 
low-albedo surfaces of the degree of redness exhibited by Hektor. By 
analogy, KBOs and Centaurs that are comparable in albedo and color 
to Hektor can be modeled with minerals (see below). 

In the models of the reddest KBOs and Centaurs described below, 
the principal component imparting the red color (increasing absorption 
toward the violet end of the reflectance spectrum) is tholin. Tholin is 
the name given to a class of complex organic solids produced by the 
irradiation of mixtures of cosmically abundant reducing gases and ice, 
and is independent of any assumptions of the detailed chemistry of 
the material. The literature of tholin production, chemistry, optical 
properties, and applicability to natural settings in the Solar System 
is slowly expanding with the growing recognition of the widespread 
occurrence of complex organic solid material in planetary environments 
(e.g., Cruikshank and Khare 2000). 

Several kinds of tholin have been made and analyzed. The starting 
mixtures, both in the form of gases and solids, have been the prin- 
cipal variables, while the energy source has typically been a plasma 
discharge producing ultraviolet light and charged particles. Table I 
shows a summary of four tholins produced at Cornell University, having 
demonstrated relevance to planetary environments. Other laboratories 
are now producing tholin materials (e.g., Ramirez et al. 2002; Imanaka 
etal. 2003). 

The Cornell University tholins, which are macromolecular complexes 
of polymers, polycyclic aromatic hydrocarbons, nitriles and nitrogen 
heterocyclic compounds (when N is included in the starting mixture), 
should be regarded as an intermediate product of the irradiation of 
simple ices and gases. As Thompson et al. (1987) showed, when they 
irradiated a methane clathrate (frozen CH 4 4- H 2 O), the initial redden- 
ing of the material eventually progressed to an overall darkening as the 
carbon was more completely dehydrogenated. Andronico etal. (1987) 
found similar results with organics synthesized in the laboratory by ion 
beam bombardment; the materials continued to darken with increas- 
ing ion fluence. Continued energy deposition in tholins is expected to 
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Table I. Tholins For Which Complex Refractive Indices Are Measured 


Name 

Starting Mixture 

Energy Source 

References 

Titan tholin 

Gaseous N 2 : CHi(9 : 1) 

Plasma Discharge 

Khareetal. 1984, 




McDonald etal. 1994 

Triton tholin 

Gaseous N 2 : CH*(99.9 : 0.1) 

Plasma Discharge 

McDonald etal. 1994 

Ice tholin 1 

H 2 O : C 2 H e {6 : 1) 

Plasma Discharge 

Khare etal. 1993, 




McDonald etal. 1996 

Ice tholin 2 

H 2 0 : CH 3 OH : CO 2 : C 2 H 6 

Plasma Discharge 

McDonald etal. 1996 


(80:16:3.2:0.8) 




remove hydrogen and other atoms, eventually resulting in a residue of 
mostly elemental carbon. 

The tholins listed in Table I are important to the present work 
because their complex refractive indices (n and k) have been measured 
over a wavelength range comparable to that in which Solar System 
bodies have been observed. Rigorous scattering models require n and k 
of the modeled constituents for the radiative transfer calculations that 
yield synthetic spectra for comparison with observational data. Very 
few materials have been measured for n and k in the spectral region 
(and with the spectral resolution) that matches the observational data 
acquired from spacecraft and ground-based telescopes. Insofar as a few 
organic solids, the tholins, have been measured and chemically char- 
acterized, these materials have proven to be a suitable and plausible 
component of the successful models calculated to date. This success 
may be fortuitous, but probably it is not. The mechanisms for energy 
deposition in simple ices and gases are prevalent in planetary environ- 
ments, and the existence of the requisite raw materials (CHON) is well 
established. 

3.2. Minerals 

In addition to ices and organic materials, some significant mineral com- 
ponent is expected in KBOs, particularly because olivine and pyroxene 
are seen in the dust of active comets (e.g., Wooden 2000). Cruikshank 
etal. (1998) found that Mg-rich olivine (Fo 82) was an important 
component in their models of 5145 Pholus, where its absorption near 
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1 jim had a significant effect in shaping the steep red continuum of this 
Centaur. 

Mafic minerals in the inner Solar System are ’’weathered” by the 
incident solar wind (e.g., Pieters etal. 2000), which has the effect of 
making their spectra substantially redder than the unaffected minerals 
exhibit. This alteration has been identified as solar wind sputtering, 
which selectively removes, reduces, and deposits the iron from mafic 
mineral grains as nanophase neutral iron on the undersides of adjacent 
grains (Hapke 2001). This process is now generally viewed as the cause 
of the red colors of Mercury, the Moon, and probably many of the main- 
belt asteroids. Processes of space weathering in the OSS are expected 
to arise from other phenomena because of the intrinsic weakness of 
the solar wind beyond several Astronomical Units (AU). For example, 
galactic cosmic rays deposit approximately 160 eV per molecule in the 
upper 100 gcm _2 of a planetary regolith at Pluto’s distance AU (John- 
son 1989), and recent estimates of the flux of anomalous cosmic rays are 
approximately a factor of 10 higher (Cooper etal. 2003). The effects 
on mineral grains in the surface zone of interaction (which may be 
tens of cm deep) can include sputtering, as well as grain amorphization 
(Brucato etal. 2003). The details of these processes in the OSS, and 
their effect (s) on the reflectance spectra of the materials in the optical 
surface remain to be studied. 

In our models we test the hypothesis that metal-reddened mafic 
minerals contribute to the red colors of Centaurs and KBOs by using 
the space-weathered lunar soil as one of the components. Shkuratov 
etal. (1999, Fig. 5) derived the optical constants for the reddened 
lunar soil by inverting their scattering model, and Dr. Shkuratov kindly 
provided the values of n and k that we use here. 

3.3. Ices 

Although ices are expected to be a major constituent of KBOs and 
Centaurs, the only ice identified with certainty on KBOs and Centaurs 
is H 2 O, seen in the spectrum through its broad absorption bands at 1.55 
and 2.0 /un , and the short wavelength part of the 2.7 /im band (e.g., 
Brown etal. 1998, 1999; Cruikshank et al. 1998; Pouletetal. 2002). 
The spectrum of 5145 Pholus shows a band complex at 2.17 /im that 
has been attributed to CH 3 OH ice (Cruikshank et al. 1998), but could 
be a light hydrocarbon. Crystalline (hexagonal) H 2 O ice can be distin- 
guished from amorphous ice, which is stable only at low temperature 
(T < 90 K) by its spectral signature, in particular the presence of a 
temperature-sensitive absorption band at 1.65 /im that is unique to 
hexagonal ice (Grundy and Schmitt, 1998). When the H 2 O ice bands 
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are only weakly present, as in the cases of the KB Os on which ice has 
been identified, it is often impossible to discern the 1.65 pm band or 
the detailed profile of the 2.0 pm band, and it is therefore unclear if 
the ice is crystalline or amorphous. 


4. MODELS OF KUIPER BELT OBJECTS AND 

CENTAURS 

4.1. Models with organics 

We have used the Shkuratov radiative transfer code (Shkuratov et ah 1999) 
for scattering in particulate surfaces to calculate models of groups II, 

III, and IV defined by Barucci etal. (2001), as shown in Figure 2. For 
each group we show three models with different proportions of the four 
basic components Titan Tholin, Ice Tholin II, amorphous carbon, and 
lunar soil. While the Barucci color data extend only to the J wavelength 
(1.22pm ), we have calculated the models to 2.5pm because additional 
color and spectroscopic data exist for some of the brighter KB Os and 
Centaurs. 

For models without the lunar soil component, we find that we can 
match the spectral gradient and the albedo level for all three groups by 
varying the proportions and grain sizes of the two tholins. The addition 
of lunar soil still allows adequate fits to the data points but introduces 
curvature to the calculated spectrum, notably a broad absorption at 
1 pm where pyroxene and feldspar in the lunar soil have overlapping 
absorption bands. The second pyroxene band around 2pm also appears 
in calculated spectra with larger amounts of lunar soil. In the reddest 
spectra, Group IV, the models allow a substantial fraction of these 
space-weathered, reddened minerals, but the absorption bands intro- 
duced into the calculated spectra are not seen in any KBO or Centaur 
spectra published to date. The model spectrum with only the tholins 
and carbon for Group IV provide a smoother fit to the data points, and 
on to 2.5pm . The situation is similar for Group III; reddened minerals 
are allowed, but produce a spectrum with broad features that are not 
seen in other data for KBOs and Centaurs. In the Group II spectra the 
absorption band structure shows much less contrast, as the calculated 
spectra begin to flatten to match the less colored objects defining this 
group. Group II admits a larger fraction of lunar soil than the previous 
groups, although the downward turn of the spectral reflectance at the 
shortest wavelengths appears to be inconsistent with the data for KBOs 
and Centaurs. 

We conclude from these models that the full range of redness of KBO 
and Centaur reflectance spectra, as well as their generally low albedos, 
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can be matched with synthetic spectra calculated from combinations 
of tholins and amorphous carbon. We find that the model spectra can 
admit, but do not require the addition of space-weathered, igneous 
rock-forming minerals. 

We noted above the photometric study of several KB Os by McBride 
et al. (2003) in which the critical V and J wavelengths were measured 
simultaneously. In Figure 3 we show our best-fit models for three repre- 
sentative objects from McBride’s list. These are models in which tholins 
provide the color and carbon adjusts the albedo. In the case of KBO 
1993 SC (top panel) we add H and K photometry from Hainaut and 
Desanti (2002) to show that the model fits the data satisfactorily. 



Figure S. Shkuratov spatial (areal) mixture and intimate mixture models of groups 
II, HI, and IV defined by Barucci etal. (2001). Key: IT2 = Ice Tholin 2, TT = 
Titan tholin, LS = space- weathered iunar soil, AC = amorphous carbon. Grain sizes 
given in pm in parentheses. Top panel: Model A is [0.36 IT2 (10) + 0.28 TT (10) 
+ 0.36 AC (10)]. Model B is [0.48 IT2 (20) + 0.32 TT (30) + 0.2 LS (20)]. Model 
C is [0.35 IT2 (20) + 0.30 TT (10) + 0.15 LS (20) + 0.20 AC (10)]. Middle panel: 
Model A is [0.25 IT2 (10) + 0.20 TT (20) + 0.55 AC (10)]. Model B is 0.28 IT2 (20) 
= 0.42 TT (60) 4- 0.30 LS (20)]. Model C is 0.28 IT2 (20) + 0.24 TT (30) + 0.12 
LS (20) + 0.36 AC (30)]. Bottom panel: Model A is 0.20 IT2 (10) ■+ 0.15 TT (30) + 
0.65 AC (10)]. Model B is [0.20 IT2 (30) + 0.06 TT (80) + 0.74 LS (30)]. Model C 
is [0.175 IT2 (30) + 0.065 TT (30) + 0.175 LS (30) + 0.585 AC (30)]. Each model 
is a spatial mix of TT and IT2, each intimately mixed with AC. 1 
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Figure 3. Shkuratov models of three objects from the McBride list. Each model is a 
spatial (area) mixture of two intimately mixed components; all of the particles are 
IQum in size, and the numbers outside the brackets indicate the fraction of the area 
covered by a specific component. Key: H 2 O = water ice at temperature 40 K, IT2 
= Ice Tholin 2, TT = Titan tholin, AC = amorphous carbon. 1993 SC: 0.45 [60% 
IT2 + 40% AC] + 0.55 [70% TT + 30% AC], 1998 SM165: 0.60 [50% IT2 + 50% 
AC] + 0.40 [70% TT + 30% AC]. 1997 CS29: 0.50 [50% IT2 + 50% AC] + 0.50 
[40% TT + 60% AC] 


Models of the kind described above can also incorporate H 2 O ice 
in amounts sufficient to display the 1.55, 2.0, and 2.5 /z m absorption 
bands. The discovery of H 2 O ice on planetary satellites, KBOs, and 
Centaurs of geometric albedo (at 0.55 /zm ) 0.1 and less (e.g., Brown 
etal. 1998, 1999; Owen etal. 1999) demonstrates that the absorp- 
tion bands are not entirely quenched by the coexistence of low-albedo 
materials on their surfaces. Indeed, in our models we find that we can 
include substantial amounts of H 2 O ice containing tholins as inclusions. 
The incorporation of tholins in the ice as inclusions is approximately 
equivalent to a mixture at the molecular level, and the optical constants 
of the mixture are calculated by effective medium theory (Shkuratov 
et al. 1999). 

Figure 4 shows two sets of models of KBO 20000 Varuna, one with- 
out H 2 O ice, and one including H 2 O ice; the actual presence of ice on 
Varuna is uncertain. The observational data for Varuna in Fig. 4 are 
from McBride et a 1. ( 2003 ). 
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Figure 4 ■ Representative models of Kuiper Bfelt object 20000 Varuna, showing a 
model without H 2 O ice (upper panel) and models including H 2 O ice in different 
abundances (lower panel). Each model is a spatial (area) mixture of two intimately 
mixed components. Key: H 2 0 = water ice at temperature 40 K, IT2 = Ice Tholin 
2, TT = Titan tholin, AC = amorphous carbon. Grain sizes are given in pm in 
parentheses, and the numbers outside the brackets indicate the fraction of the area 
covered by a specific component. Upper panel: spatial mix of 60% of component 1, 
which is an intimate mix of [50% IT2 (10) + 50% AC (10)], and 40% of component 
2, which is an intimate mix of [60% TT (10) 4- 0.4 AC (10)]. Lower panel: Model 
A is 0.5 [30% H 2 O (10) with inclusions amounting to 50% IT2, intimately mixed 
with 70% AC (15)] + 0.5 [80% H 2 O (10) with inclusions amounting to 50% TT, 
intimately mixed with 20% AC (15)]. Model B is 0.65 [50% H 2 0 (30) with inclusions 
amounting to 50% IT2, intimately mixed with 50% AC (30)] + 0.35 [50% H 2 O (30) 
with inclusions amounting to 50% TT, intimately mixed with 50% AC (30)] 

4.2. Models without organics 

Objects having a moderate degree of redness, such as the D-type as- 
teroids (mostly Trojans asteroids) and the objects in Groups I and II 
of Barucci etal. (2001) do not require the presence of organic mate- 
rials to achieve their red colors. Cruikshank et al. (2002) modeled the 
spectrum of Trojan asteroid 624 Hektor with simple Hapke scattering 
models using only mafic minerals (Mg-rich pyroxene) and amorphous 
carbon. Their models included the spectral region out to A = 3.6 /im , 
which covers the region of absorption bands in hydrous silicates; The 
models demonstrated that a significant amount of hydrous silicates 
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(e.g., serpentine) can be present on Hektor without causing a ab- 
sorption band that is detectable with data of modest quality. Figure 
5, adapted from Cruikshank etal. (2001) shows the data for Hektor 
and a suite of spectral reflectance models using intimate mixtures of 
pyroxene, amorphous carbon, and serpentine. As before, the low albedo 
is achieved by the use of the carbon, while in these models the pyroxene 
imparts the red color. In the models in the bottom panel of Fig. 4, 
graphite is included; this component has a slightly blue color at the 
long wavelengths. 



Figure 5. Models of D-type asteroid 624 Hektor without organics. Pyx80 = pyrox- 
ene with 0.8 Mg, Pyx70 = pyroxene with 0.7 Mg, AC = amorphous carbon, Serp 
= serpentine, Graph = graphite. Grain sizes are given in parentheses. These are 
all intimate mixtures, with the proportions shown as percentages. Adapted from 
Cruikshank et al. (2001) 
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5. CONCLUSIONS 

We find that we can match the geometric albedos and colors of the 
reddest color groupings of outer Solar System bodies with models com- 
puted with scattering theory and the use of optical properties of syn- 
thetic organic solids (tholins), particularly Ice Tholin II and Titan 
tholin. In our models, elemental amorphous carbon must be added 
to the organic solids to achieve the low geometric albedos observed 
or assumed for these objects. The presence of elemental carbon can 
occur in nature by the dehydrogenation of tholins. In the use of tholins 
to achieve the observed red colors, our results are consistent with 
those of several other investigators, including the recent work of Dotto 
etal. (2003) and Doressoundiram etal. (2003). The least red objects 
in the outer Solar System can be modeled with mafic minerals, and do 
not require the presence of organic solids (Cruikshank etal. 2001). 

We have incorporated into our models the space- weathered minerals 
found on the lunar surface. The Moon, Mercury, and many asteroids 
are reddened by sputtering of iron from iron-bearing minerals, and 
the deposition of neutral Fe on grains in the uppermost regolith. Our 
models of the reddest bodies in the outer Solar System do not require 
this mechanism to achieve the observed albedos and colors, although 
such space-weathered minerals may occur in modest quantities in the 
optically accessible surfaces. 
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Notes 

1 The number of significant figures in some of the modeling parameters suggests 
a degree of sensitivity of the fits to the data that is perhaps misleading. While the 
quality of fit of the preferred models is in fact sensitive to these parameters at a high 
level of accuracy, the models themselves have many more physical parameters that 
are not well known for most bodies in the Solar System. In addition, the complex 
refractive indices of the materials we use in the models are in some cases uncertain 
by 30% or more. Once a set of refractive indices is selected and other parameters 
of particle scattering are chosen, the calculated models are quite sensitive to the 
abundances of the components, their grain sizes, and their mixing geometry. 
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